Complex Acquisition of the Fourier Transform Imaging of an Arbitrary Object 
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A scheme to a complex-valued acquisition of the Fourier transform imaging was proposed. The 
main idea is to project the real and the imaginary parts of a diffraction field to intensity distributions 
respectively. The whole procedure was algorithm independent and needs no a priori knowledge of 
an arbitrary objet. An example was demonstrated with a numerical modeling and its results. 
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When talk about the image forming, the process by 
instruments like the eye, the camera, the reflecting and 
refracting telescope, and the microscope etc., were gen- 
erally refers to a point to point correspondence between 
two real spaces |1] . If a situation requires full information 
about the objects rather than only a magnitude transmit- 
tance or reflectance, such procedures would be invalid be- 
cause what they mapped was only intensity relations, and 
phase information would thus lost. A possible way to ob- 
tain the phase knowledge of the object's transmittance 
is to convert it into a spatially varying diffraction pat- 
tern with form of Fourier transform (supposed in Fraun- 
hofcr region), and then, invert it to its object function. 
This procedure is based on the fact that Fourier trans- 
form keeps unitary relations between real and recipro- 
cal spaces. Unfortunately, the diffraction fields are also 
only recordable by an intensity-sensitive detector and this 
awkward fact would lead phase loss to occur in the re- 
ciprocal space again. To solve the phase problem, ef- 
forts have been paid by utilizingoversampling methods 
with iterative algorithms [^,111, 0]. The modulus of 
Fourier transform can thus be phased and then inverting 
it into an object functions. This method was reported 
recently in X-ray diffraction microscopy ^ and has been 
extended from x-ray crystallography 7] to the imaging of 
noncrystalline materials [1], and single cells [9]. Its po- 
tential for imaging of single protein complexes by using 
ultra-short X-ray pulses with extreme intensity were also 
discussed [l^l- The excellent works as mentioned above 
have now achieved the imaging of single virion even with 
a resolution of 22nm (Tlj . 

Anyway, the phase-less magnitude alone would not be 
able to sustain the unitary property of Fourier transform. 
Therefore, iterative algorithms to recover the phase en- 
coded in the diffraction pattern have to rely on a pri- 
ori knowledge of the objects more or less. If a object 
was arbitrary i.e. not purely absomtive or not purely 
phased, the ambiguity would arise [ij]- The difficulties 
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FIG. 1: The proposed scheme for set up 



of the complex-valued acquisitions without any original 
information seemed to have already been predestined by 
quantum mechanics, i.e., the complex amplitude can not 
be specified exactly in a single measurement |13i] . On the 
other hand, for coherent sources are not obtainable in 
wave length like hard X-ray region, potential advantages 
might be jeopardized when using the ultra-short, in- 
tense X-ray pulses jTo|. 

In this paper, we theoretically propose a methods and 
a scheme to a complex- valued acquisition of the Fourier 
transform imaging of the arbitrary object's transmittance 
with incoherent light. The main thought is to project the 
real and the imaginary parts of a Fourier-transformed 
field to intensity distributions respectively. The whole 
procedure is algorithms free and needs no a priori knowl- 
edge of an arbitrary objet. 

The proposed scheme for set up is shown in Fig[T] 
Fields from the thermal source is split by a 50/50 beam 
splitter BSi to form a two-arm optical system with equal 
distance from source plane x to planes rji, and r]2, where 
intensity information will be registered. At the cross sec- 
tion of both arms before rji, and 772, the other 50/50 beam 
splitter BS2 would insert. As for BS2, we must manage 
to ensure that the half wave loss only occurs when the 
light reflects only on one of its surfaces. When a situation 
needs to introduce a phase factor shift of j, a phase plate 
J would be inserted into the upper part of the system. 
The total scheme differs mainly from the Mach-Zehnder 
interferometer is that it uses a prismatic lens P rather 
than a plane mirror to guide the optical path in its up- 
per part. The object with complex transmittance of /(^) 
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is placed at plane ^ in lower part of the setup. The dis- 
tance from plane ^ to plane x, and to plane 771, and from 
plane x to plane 772 are c?i, ^2 and d respectively. The 
equal length of the two arms requires that 



d — di+ d2 



(1) 



Among the setup, an optional phase plate P' might insert 
at the lower part of the scheme to form a fixed optical 
path difference between two arms of a phase factor of 
l/v7; whereas, it was not a necessary component. 

Although chaotic light fields fluctuates randomly, the 
relations among instantaneous values of optical fields E 
in plane ryi , 772 and ^ are deterministic according to Fres- 
nel diffraction theory. To illustrate the theoretical bases 
for this complex-valued retrieval procedure clearly, first 
we suppose, that we did not facilitate the beam splitter 
BS2, phase plate P', and J in FiglTJ In this situation, 
the Optical fields on 771 , and 771 fulfill [Tsj 
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E{m) = ^ I E{x)e^''^^dx, 
jXd 
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when half wave loss only on one side of BS2 was con- 
sidered. Their corresponding intensity distribution reg- 
istered on plane 771, and 772 are 



and 



him) -{El im)E*Avi)) = l{\EMf 

-2Re {E{m)E* (772)) + 1^(772) 



i2iv2)-{E2im)E*Am)) 



(8) 



1{\EM\' 



+2Re {E{m)E* (772)) + 1^(772)1- 



(9) 



respectively. The fact can be easily seen from two equa- 
tions above that the real part of Eq. ^ is embedded in the 
intensity registration on both plane 7^1, and 772, and can 
be extracted by subtracting the two equal backgrounds 
by means of 
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and 
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under the Fresnel approximation. Let the thermal source 
is totally chaotic with evenly distributed intensity, 
i.e., 



{E{x)E{x')) ^15{x-x'), 
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Further more, the linearity Eq.(l5]) indicates that if a 
phase shift of (/> is introduced in the upper or the lower 
part of the scheme shown as FiglU a phase factor of e~^'^ 
or e^'^ would be multiplied on the right side of the equa- 
tion consequently. Following the thought, if we inserted 
a phase plate J into the upper part of the system to in- 
troduce a phase shift of 7r/2, the complexes vector stands 
for Eq.® would rotate an angle of —tt/2 consequently. 
So the imaginary part of Eq. ^ can also be retrieved in 
a way similar to eg. pop by: 
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in which (• • •) stands for the assemble average. Then, the 
derivation of the mutual intensity between plane 771 and 
ri2 can be derived from Eq. H]) to Eq.® to be 



{E{m)E* (7,2)) = v^e 



XF I — 1 CX y/je 



I 



\d2 



■K\J X^ddidl 

^{fim- (5) 



In which F!F{- ■ ■) refers to JF{/(a;)},thc Fourier transform 
of f{x). Now with the beam splitter BS2 facilitated, the 
fields in plane 771, and 772 turn out to be: 



and 
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since a factor of —j had been brought into Eq.®. In 
Eq. pT]) . 12(112) and /i(77i) are intensity registration on 
plane 772 , and 771 , after a phase plate J was inserted into 
the upper part of the system as FiglT] shows. Note that 
we use a prismatic lens P rather than a plane mirror 
to guide the optical path in upper part of the system. 
This arrangement leads to a bilateral symmetry between 
coordinates in plane 771, and 772, i.e., 



vi = -m (= ri) ■ 



(12) 



Comparing Eq.®, Eq.dTO]), Eq.lTTl), and Eq.((T2l), we 
proposed the complex-valued acquisition of an object's 
Fourier transform imaging. The procedure can be writ- 
ten in one equation as: 
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FIG. 2: The real (a) and imaginary (b) parts of the Fourier 
transform as a function of the space coordinates rj. It origi- 
nates from a complex- valued transmittance (eq.(14)) of a con- 
ceived object. 



To give an example of the retrieval, we conceived an 
object with a complex-valued transmittance of: 
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In which 0.05, 150, 105, and 1000 are space parame- 
ters with unit in plane ^. If the object were to 
be illuminated by a coherent light with a wavelength of 
A — 0.532/^TO, the real and imaginary parts of its Fourier 
transform as a function of space coordinates rj would be: 



Re (J^ {/(C)}) = 500 sine 
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and 



lva{T{J (e)}) = 210sine(l05^j 
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as FigHihows. In Eq.dini) and Eq.lHD, = 75,000^m, 
refers to a distance the coherent light field with A = 
0.532/im from the object propagates until it reaches the 
plane rj. 

Based on the theory of statistical optics, we numeri- 
cally modeled the dynamic process of the whole retrieval 
procedure under the setup scheme of FiglU by using the 
conceived object with transmittance of Eq. (fT4|) . An op- 
tional phase plate prism P' was inserted at the lower part 
of the scheme to form a fixed optical path difference of 
7r/4 between two arms to compensate the factor of 1/a/J 
in Eq.®. 

In the modeling, the monochromatic thermal light fea- 
tured circular Gaussian random process with zero mean 
, the wave length of the thermal light was selected to 
be A = 0.532/xm, the propagation of the fields among the 
setup fulfills the Fresnel diffraction integrations stated 




FIG. 3: Numerical results for the Complex- valued Retrieval 
of a conceived object. Its left and right column concerns in- 
tensity registration before and after a phase plate J of 7r/2 
was inserted in the upper part of scheme in Fig[T] Comparing 
(c) and (c') with (a) and (b) in FIG.[2l the numerical results 
shows that the real and imaginary parts of the Fourier imag- 
ing of a conceived object (Ea. (|14p l were both retrieved in a 
sub-wave length scale. 



by Eq.(I2)to Eq.®. For the setup, di, ^2 and d are 
set to 60, OOO/Ltm, 75, 000/xm, and 135, 000/zm. After an 
accumulative intensity registration, which covers 20, 000 
times of independent coherent time, the averaged inten- 
sities in both plane rji and 772 were brought about in FIG. 
3. In which (a) and (a') are intensity registered in plane 
ri2 before and after a phase plate of 7r/2, i.e., J was in- 
serted; (b) and (b') are intensity registered in plane rji 
before and after a phase plate of J was inserted; (c) and 
(c') was the intensity difference between plane 772 and 771 
before and after J was inserted. Comparing (c) and (c') 
with (a) and (b) in Figl^l the numerical results shows 
that the real and imaginary parts of the Fourier imaging 
of a conceived object fEg. pi)) ) were both retrieved in a 
double scale coordinate of FigO 

In summary, by using incoherent light, we theoreti- 
cally proposed a scheme to a complex acquisition of the 
Fourier transform imaging. The whole procedure needs 
no a priori knowledge of an arbitrary object. The unitary 
property of the Fourier transform thus sustains by inten- 
sity recordable detectors. In recent years, coincidence 
imaging by classical thermal sources has been widely 
investigated both in real [l3, [3 ^'^^ reciprocal spaces 
[T9I [20 . [2H . In these works, images were carried out 
either by two photon coincident rate or by correlation 
functions of intensity fluctuations. They both stand for 
the modules of optical fields. Useful methods to retrieval 
a complex-valued objects with discrete variables [l^, [2^ 
were also reported combined with oversampling meth- 
ods [1| and iterative algorithms [1, 0, The authors 
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of Ref.fl^l reconstruct the complex object without a pri- 
ori knowledge, but still based on computer operations. 
Unlike these works, we map the complex valued object 
in real space to complex valued Fourier transform only 
by intensity registration (23 |. As the unitary property 
remains, the retrieval procedure by this scheme equals 



to imaging an arbitrary complex valued object in a full 
sense. To our knowledge, it is the first physical proposal 
to sustain the unitary property of an arbitrary object 
when imaging without a priori knowledge. 
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